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We consider second-order elliptic partial differential operators acting on sec- 
tions of vector bundles over a compact Riemannian manifold without bound- 
ary, working without the assumption of Laplace-like principal part — V^V^ . 
Our objective is to obtain information on the asymptotic expansions of the 
corresponding resolvent and the heat kernel. The heat kernel and the Green's 
function are constructed explicitly in the leading order. The first two coef- 
ficients of the heat kernel asymptotic expansion are computed explicitly. A 
new semi-classical ansatz as well as the complete recursion system for the heat 
kernel of non-Laplace type operators is constructed. Some particular cases are 
studied in more detail. 



Ivan Avramidi and Thomas Branson: Non-Laplace type operators 



1 



1 Introduction 

The resolvent and the heat kernel of elliptic partial differential operators are of great 



importance in mathematical physics, differential geometry and quantum theory |L7, [L2 



If, 19 1 . Of special interest in the study of elliptic operators are the so-called heat kernel 



asymptotics. It is well known [12| that for a second-order, elliptic, self-adjoint partial 
differential operator F with a positive definite leading symbol, acting on sections of a 
vector bundle over a compact, boundariless manifold M of dimension m, an asymptotic 
expansion of the following form is valid as t I 0: 



Tr L2 exp(-tF) ~ (A-Kt)'^ 2 £ -^-J—A k . 



k>0 



The coefficients A k are called the heat invariants, or heat kernel coefficients. They are 
spectral invariants of the operator F, and encode information about the asymptotic prop- 
erties of the spectrum. Note that our normalization of the coefficients A k differs by the 
factor (— l) k /k\ from that used in [12]. This normalization has been used in previous 
works of one of the authors (see the book ||, the papers @, §, §] and others). It has 
the advantage that for Laplace type operators with a potential (see definition below), 
F = —A + q, the numerical coefficient of the term j M dwo\ (x) tr q k in is equal to 1 for 
any k. 

An important subclass of the class of operators described above are the operators 
of Laplace type: those for which the leading symbol takes the form g^C.^C.u , where g^ u 
is a non-degenerate, positive definite metric on the cotangent bundle of M. For such 
operators, the leading symbol naturally defines a Riemannian metric on M (the inverse g^ v 
of the leading symbol). Alternatively, one may take the Riemannian metric as given, and 
produce the many natural operators of Laplace type which are so important in Physics. 

The assumption of Laplace type affords a considerable simplification in the study 
of spectral asymptotics. Partly as a result of this, much is known about the resolvent, 
the heat kernel, and the zeta function in this category of operators |12|, || [|, |], In 
particular, the heat kernel coefficients for Laplace type operators are known explicitly 
up to k = 4 @, H (for a review, see Q). 

In this paper, we take a Riemannian metric as given, and study the most general class 
of second-order operators F, acting on sections of a vector bundle V, with positive definite 
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leading symbol. That is, we drop the assumption of Laplace type, and assume only that 

a 2 (F;x,0 = a^(x)^ v , 

where a? v is a symmetric two-tensor valued in End(V). (We do not assume that = 
g^idy , nor that a^ u is factored as g^ u E for E a section of Aut(V).) We shall sometimes 
call these "NLT" (for "non-Laplace type" ) operators. Of course, Laplace type is a special 
case, so to be more precise, NLT operators are operators that are not necessarily of Laplace 
type. NLT operators arise naturally in such areas of mathematical physics as quantum 
gauge field theory and quantum gravity [[II], |J, differential geometry, classical continuum 



mechanics |18| and others. 

The most elementary examples that one can use to illustrate the class of operators in 
question are the weighted form Laplacians. If d is the exterior derivative and 5 its formal 
adjoint, the form Laplacian is A = Sd + dS; this differs from the Bochner Laplacian 
fl^V^Vj, by an operator of order zero, the so-called Bochner- Weitzenbock operator. If 
a, b are real constants, an operator of the form aSd + bd5 may be termed a weighted form 
Laplacian; such operators are elliptic if a ^ ^ b, and have positive definite leading 
symbol if a, b > 0. 

The non-Laplace type operators on differential forms have been studied extensively 
by many authors under the general name "nonminimal operators" , or "exotic operators" 
(see |T3|, []. 113, [T3|, [Tj], [TDL |T[] and references therein). In particular, the paper |T3|] contains 
a complete discussion of the coefficients of the weighted form Laplacian for all k and of 
A and Ai of the operator with a potential function (Th. 1.2 and Th 1.3, pp. 2089-2090). 
In papers fl4], [l^, [16[ a computer algorithm has been developed that employs the calculus 



of pseudodifferential operators and the coefficients A , Ai and A 2 (more precisely, the 
local (non-itegrated) coefficients a , a\ and a 2 ) have been computed. 

The ^-function for non-Laplace type operators acting on symmetric 2-tensors have been 
studied in |I0| (for the origin of such operators, see [11], |5|). By restricting to maximally 



symmetric even dimensional manifolds the authors were able to compute the eigenvalues 
of such operators and evaluate the ^-function at zero, C(0)> (which is, in fact, equivalent 
to the coefficient A m / 2 fL2|| ) for a certain special case of the operator in Euclidean spaces 
and m-spheres in dimensions m = 2,4, 6, 8, 10. 

In general, the study of spin-tensor quantum gauge fields in arbitrary gauge necessarily 
leads to non-Laplace type operators acting on sections of general tensor-spinor bundles. 



Ivan Avramidi and Thomas Branson: Non-Laplace type operators 



3 



It is precisely these operators that are of prime interest in the present paper. The main 
examples here are the symmetric 2-tensor bundle (gravitational field, spin 2) and the 
spin- vector bundle (gravitino field, spin 3/2) (for a discussion of such operators in gauge 
field theories, see |7|]). 

Let us formulate our main result from the very beginning. 

Theorem 1 Let F : C°°(V) — > C°°(V) be a self-adjoint elliptic second-order partial 
differential operator with a positive definite leading symbol, acting on sections of a tensor- 
spinor bundle V of fiber dimension d over a compact manifold M of dimension m without 
boundary. Let F = — V M V^ + q, where V is a connection on the vector bundle V , a^ v 
is a parallel symmetric two-tensor valued in End (V) and q is an endomorphism of the 
bundle V. The curvature of the connection V is defined by [V M , V u ]<p = R a p^T^ a (p, where 
R a piiv is the Riemann curvature tensor and T a p is given by the representation of so(m) 
which induces the bundle V. Let £ G I 'M be a cotangent vector, and let Aj(£) = /ij|£| 2 , 
(i = l,---,s), fii > 0, be the eigenvalues of the leading symbol, (72(F) = a^^^u, with 
the multiplicities d{. Then the corresponding (orthogonal) eigenspace projections have the 
form 



k=i 

where the numbers s and p depend on the structure of the bundle V and the leading symbol. 
The matrix of coefficients C := (cik) is inverse to the (Vandermonde) matrix of powers 



Furthermore, the L 2 trace of the heat kernel has the following asymptotics as t — > 




n=0 




M = (n kj ) := Ou*- 1 ). 



Tr L 2 exp(-tF) = (Aut)- m ' 2 [A - tA 1 + 0(t 2 )] . 



The coefficients Aq and A\ are defined by 



^o = ]T^/i- m/2 vol(M), 



A 



1 — 



d vol (x) { tr v (a Q q) + (3 R} , 



M 
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where R is the scalar curvature, 



s 



i=i 



^ n ^ n r(m/2)(2fc-2)! 



k-- 



^ r(m/2 + k — l)2 2k ~ 2 (k — 1)! 



and (3 is a constant defined by 



p = -^E^ rm/2tr v(<^>^o 



i=i 



+ 10 , 1 V [Kifc(3m - 2) + 4(m + 2)^<r ifc ] tr v (rL;J a n fe J Q 

12(m — 1) , \ 

' l<i,k<s;i^k 



with 



2(m- 1) 

l<i,fc<s;i^fc 



r(m/2 - 1) f /i- 72 + /^r /2 , o ^~ m/2+1 ~ h m ' 2+1 
lk 2r(m/2 + 1) 1 - " - /i fc ) 2 



-m/2 . -m/2 

8r(m/2 + 2)i" (/i.-Zifc) 3 v (^-^ fc ) 2 

-m/2 -m/2-1 

+4(m-2)— ^ - + m(m-2) ^ 



(/ii-/ifc) 2 ' {Hi- Hk) y 



j j sr r(m/2)(2n + 2j-2)! 

^ > i 2^ r(m/2 + n + j - l)2^-\n + j - 1)! 

^9(fllfJ,2 ' ' ' 9fJ.2n-[ifl.2n-29vil/2 ' ' ' 9u 2 j-^2j-29^/S) 
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We shall prove this theorem in Section |9|. Note that for Laplace type operators, when 
a^ u = g^Iy, these formulas simplify considerably. We have then just one eigenvalue 
Hx — 1 with multiplicity d\ = d = dim(V) and the projection III = Iy. Thus for Laplace 
type operators, ao =< IT >= iy, and /3 = —d/6, and we recover the well known result 

A = vol(M), 
A x = J dvo\(x) ^tr v g- 

In the next section we shall give a detailed description of our operator class. We would 
like to stress from the beginning that the theory of NLT operators, despite being a theory 
of second order operators, is closely related to the theory of higher-order Laplace-like 
operators (for a discussion, see ||); that is, operators with principal part a power of the 
Laplacian. In this sense, even the weighted form Laplacian example mentioned above 
does not capture the full flavor of the theory, since it (almost uniquely in this category) 
is accessible entirely through second-order methods (see |13|, jlBf). 

To contrast with the previous result mentioned above, it is worth stressing once again 
that: 

i) there are many areas in both physics and mathematics where general second-order 
non-Laplace type operators arise naturally, 

ii) we are primarily interested not just in the weighted form Laplacians but in the 
general NLT operators (differential forms being a very particular case), 

iii) our approach to computation of heat kernel asymptotics is completely different from 



that of the previous authors |13|, |T6|| . Our method enables us to compute explicitly 
not just the heat trace asymptotics but also the heat kernel and resolvent in the 
leading order that describe the local off-diagonal behavior of the resolvent and the 
heat kernel. To best of our knowledge, such formulas for non-Laplace operators are 
presented here for the first time. 

Despite the importance of second-order operators with non-Laplace principal part in 
gauge field theory and quantum gravity (see, for example |Ll], ||), their study is still 
quite new, and the available methodology is still underdeveloped in comparison with the 
Laplace type theory. In this paper, and in the more explicitly representation-theoretic 
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treatment 0, we hope to lay the groundwork for a systematic attack on the spectral 
asymptotics of this larger class of operators. 

2 Non-Laplace Type Differential Operators 
2.1 General vector bundle setup 

Let M be a smooth compact manifold without boundary of dimension m, equipped with 
a (positive definite) Riemannian metric g. Let V be a smooth vector bundle over M, with 
End(V) = V <S> V* the corresponding bundle of endomorphisms. Given any vector bundle 
V, we denote by C°°(M, V), or just C°°(V), its space of smooth sections. We assume that 
the vector bundle V is equipped with a Hermitian metric H. This naturally identifies the 
dual vector bundle V* with V, and defines a natural L 2 inner product, using the invariant 
Riemannian measure dvol(x) on the manifold M. The completion of C°°(M, V) in this 
norm defines the Hilbert space L 2 (M, V) of square integrable sections. 

We denote by TM and T*M the tangent and cotangent bundles of M. We assume 
given a connection V v on the vector bundle V. The covariant derivative on V is then a 
map 

V v : C°°(y) -> C°°{T*M® V) (2.1) 

which we assume to be compatible with the Hermitian metric on the vector bundle V, 
in the sense that V-ff" = 0. Here the connection is given its unique natural extension 
to bundles in the tensor algebra over V and V*; in particular, to the bundle V* ® V* 
of which if is a section. In fact, using the Levi-Civita connection V LC of the metric g 
together with V v , we naturally obtain connections on all bundles in the tensor algebra 
over V, V*, TM, T*M; the resulting connection will usually be denoted just by V. It will 
usually be clear which bundle's connection is being referred to, from the nature of the 
section being acted upon. We denote the curvature of V v (a section of T*M <S> T*M £g> V) 
by K: 

[V a ,V p ]ip = K af3 ip, <peC°°(V). (2.2) 
The formal adjoint of the covariant derivative of ( j2.1|) is defined using the Riemannian 



Ivan Avramidi and Thomas Branson: Non-Laplace type operators 



7 



metric and the Hermitian structure on V: 

(V v )* : C°°(T*M(g) V) ^C°°(V), 

Let 

a G C°°(TM <g> TM) <g> End(V)), b G C°°(TM <g> V), g G C°°(End(V)). 

These sections define certain natural bundle maps by contraction, which, for simplicity, 
we denote by the same letters: 

a: T*M ®V ^ TM ® V, (2.3) 

6: V^TM®V, 
g: V^V. 

Using these maps we can write the general second order operator 

F = V*(aV)+6V + g, where a^ = a u ". (2.4) 
Any formally self-adjoint second order operator may thus be written 

V(aV) + ^(6V + (6Vn+?, (2.5) 

where we may assume that each a? ,v is Hermitian. (Here, if necessary, we clear the notation 
and redefine q and a^ v '.) In abstract index notation, any formally self-adjoint second order 
operator may be written 

F = -V M (a^V,) + 6^V M - V M (6*r + g 

= -a^V M V, + [&" - (6*)" - a^, v ]V, + q - (&*)% 

Hereafter we denote by " ;" the covariant derivative. One may now restrict to the case when 
the endomorphism b* is anti-Hermitian, since the Hermitian part of b only contributes at 
order zero. With this, we have redefined b and q. Thus, henceforth, 
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The formal adjoint to the operator F reads 

F* = V*(a*V) + 6V + (6V)* + q*. 

Hence, in addition to the other conditions posited so far, for the operator F to be formally 
self-adjoint the endomorphism q should be Hermitian. To sum up, the general formally 



self-adjoint second order operator is described by (|2.5|) with 



a"" = a^, (a"")* = a"", (&")* = q* = q. 

Let us consider the effect of a change of the connection, V — > V = V + A, with a one 
form A valued in End (V). We have 

F = -V^Vu + (6 M - A,a^)V M + V M (6 M - aT A v ) + b^A^ + AJf - A^ v A v + q. 

In many cases (but not always!) it is possible to choose A in such a way that A^a^ v = 
W . Then the first order part drops out. The point is whether the map a ( |2.3|) is invertible, 
i.e. whether there is a solution, cT 1 G C°°(T*M <g> T*M) ® End(V)), to the equation 

a^a-\ x = 6£I V . (2.6) 

This can be put in another form. Let e» G C°°(T*M (g) V) be the basis in the space of one 
forms valued in V and e* G C°°(T*M ® V*) be the adjoint basis in the space of one forms 
valued in V*. Then the equation (|2.6|) has a unique solution if and only if the bilinear 
form By =< e*, ati > is nondegenerate, i.e. det B^ ^ 0. If this condition is satisfied then 
one can always redefine the connection in such a way, viz. A^ = b v a~ x ViJ ,, that A^a^ = b v 
and the first order terms are not present. In this paper we assume that this is the case, 
so that without loss of generality one can set the vector-endomorphism b to zero, 6 = 0. 
Moreover, we will assume that the tensor-endomorphism a is parallel, Va = 0. Thus the 
operator under consideration has the form 

F = -arVpVv + q, (2.7) 

where 

a^ = a ^, a*"" = a"", Va = 0, q* = q. (2.8) 
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2.2 Tensor-spinor bundles 

We now restrict attention to operators acting on tensor-spinor bundles. These bundles 
may be characterized as those appearing as direct summands of iterated tensor products 
of the cotangent and spinor bundles. Alternatively, they may be described abstractly as 
bundles associated to representations of 0(m), SO(m), Spin(m), or Pin(m), depending 
on how much structure we assume of our manifold. These are extremely interesting 
and important bundles, as they describe the fields in Euclidean quantum field theory. 
More general bundles appearing in field theory are actually tensor products of these with 
auxiliary bundles, usually carrying another (gauge) group structure. The connection on 
the tensor-spinor bundles is built in a canonical way from the Levi-Civita connection and 
its curvature is: 

<t-> _ pa rp/3 
'^■fiu J ^ P^v ± at 

where is the Riemann curvature tensor, and T a p is determined by the representation 

of so(m) which induces the bundle V. T^ a is a tensor-spinor constructed purely from 
Kronecker symbols, together with the fundamental tensor-spinor 7^ if spin structure is 
involved. 

We study in this paper a special class of second-order operators of the form ( |2.7| ) with 
the coefficient a built in a universal, polynomial way, using tensor product and contraction 
from the metric g and its inverse g*, together with (if applicable) the volume form E 
and/or the fundamental tensor-spinor 7. Such a tensor-endomorphism a is obviously 
parallel. (E is available given SO(m) or Spin(m) structure; 7 is available given Spin(m) 
or Pin(m) structure.) We do not set any conditions on the endomorphism q, except that 
it should be Hermitian. An important subclass of this class of operators is the class of 
natural operators, when, in addition, q is also built from the geometric invariants only; 
i.e. from g, g*, E, and 7, together with the Riemann curvature and its iterated covariant 
derivatives. By Weyl's invariant theory and dimensional analysis (i.e., a check of the 
homogeneity of each term under uniform dilation of the metric), it is clear that while a 
must be built polynomially from g and g*, together with E and/or 7 if applicable, the 
endomorphism q must be a sum of terms linear in the curvature. However, we do not 
need the additional assumption that q is constructed from the curvature. In general, it 
could be any smooth endomorphism. 
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3 Leading Symbol of an NLT Operator 

Let us describe now more exactly the class of operators ( ]2.7|) we are working with. We 
have assumed the operator F to be self-adjoint leading to the conditions (fDf ). Since 
we are going to study the heat kernel asymptotics of the operator F, we now require in 
addition that the leading symbol of the operator F, 

<r 2 (F)(0 =: A(£) = a^ v , with £ G T*M, 

be positive definite, i.e. we have 

£ 7^ =>- A(£) Hermitian and positive definite on V. 

In particular, F is elliptic. Positive definiteness implies that the roots of the characteristic 
polynomial 

X.(fl(A):= det v (A(£)-A) (3.1) 

are positive functions on M. 

Remark 1 An important point is that the eigenvalues Ai , . . . , \ s and their multiplicities 
d\ , . . . , d s are independent of the point x G M. Indeed, all the bundles under con- 
sideration are vector bundles associated to the principal bundle S of spin frames via a 
finite-dimensional representation {ip, V) of Spin(m): V = S x v V. Let V\ be the defining 
representation of SO(ra). We get the section a, and the endomorphisms a^ u ^^ v by "pro- 
moting" vectors a G V\ ® V\ ® V <g) V * and a^^^iv G V <8> V* to sections of the associated 
bundles. In particular, the cited eigenvalues and multiplicities may be computed at the 
level of the representation (<p, V). 

A useful way in which to perturb our operator F is to let it run through a one- 
parameter family F{e) for which 

a(e) = a + eb, 

where b is a section of TM®TM®End(V) which is built from g and g*, and if applicable, E 
and/or 7. (This b is not to be confused with the b of section 2.1.) To preserve formal self- 
adjointness, we need to assume that lf v ^^ v is self-adjoint on V for each £ G C°° (T*M). 

F = -( a » u +eb^)S7 l X v + q. 
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For e in some interval about 0, the symbol a(e) remains positive definite. (This last 
statement does not even depend on the compactness of M; by an argument analogous 
to Remark [I], this interval is independent of the point x G M.) In particular, we might 
take perturbations about an operator with a leading symbol which is factored; i.e., one 
for which 

a"" = g^c, 

where c is a section of End(V) which is built invariantly from g and g*, and if applicable, 
E and/or 7. As a special case of this, we could take c = I\> ; that is, perturb an operator 
of Laplace type in NLT directions. In fact, in case V is associated to an irreducible 
representation of Spin(m), the endomorphism c must be a multiple of the identity Iy by 
Schur's Lemma. 

Given such a perturbation, one might hope that relevant spectral quantities could be 
expanded in powers of (very small) e, or at least that one could work with the e-variation 
of such quantities. 

Now consider the symmetric 2n-tensor quantity 

tr v a^ lUl ■■■a^ Un \ (3.2) 

As usual, the parentheses denote complete symmetrization over all included indices. An 
index-free way of writing this is 

tr v V n a, 

where V n is the symmetrized tensor power defined by 

V n a = a V ■ ■ • V a y . 

n 

Remark 2 We claim that 

tr v V n a = a (n) V n g* (3.3) 
with some constants ai n ) ■ (Since Va and Vg vanish, so must Vay .) 

Indeed, this quantity is built polynomially from g and g*, and a priori possibly E 
and/or 7. Since it is a tensor, 7 is not involved. Since 

E^-^E Ul .. Mm = 5^ [ui ...S^ Umh 
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where the brackets denote antisymmetrization over enclosed indices, ( |3.2j ) is affine linear 
in E; that is, it has the form 

(ii\v\ ...fi n u n ) I I (fiiui ...n„v n )[ai ...a m ] tp 
Y 1" iy J^cei ...a m - 

The tensor ip is constructed purely from the metric g* and is symmetric in the first 2n 
indices and antisymmetric in the last m indices. It is clear that there are no such tensors, 
so ijj — 0. Similarly, the only symmetric 2n-tensor constructed from the metric is V n g*, 
thus leading to (|3.3|). 

Contracting with £^ Ul . . . ^ n ^ n , we get 



tr v A"(0=a (n) |e| 



2n 



Furthermore, if we denote by tr 5 the total trace of a symmetric 2n-tensor P, 



then since 



T(m/2 + n) 2 2n n\ 



(3.4) 



tr„ V l g 



we have from ( |3.3| ) 



t(n) 



r(m/2) (2n)! 
r(m/2) (2n)! 



tr v tr 9 V n a. 



r(m/2 + n) 2 2n n\ 

It is clear that for an operator of Laplace type, when a = g* ® idy, each a( n ) is just the 
fiber dimension d of V. 

Now let us consider the characteristic polynomial (|3.1|) in more detail. Applying the 
above remarks to the coefficients of Xa(0 (A), we find that Xa(0 depends on £ only through 
|£| 2 . As a result, the dependence of the eigenvalues Aj on £ is only through |£| 2 . Since 
A(£) is 2-homogeneous in £, the Aj must be also: 

for some positive real numbers \i\ , . . . which are independent of the point G 
T*M, and, in fact, independent of the specific Riemannian manifold (M,g), depending 
only on the representation (ip, V) and the vector a e Vi <g> ® V" ® V*. Computing the 
trace trv^4(0 n we § e t a sequence of equations, 



(3.5) 
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relating the eigenvalues, multiplicities, and the quantities a( n ) . When n — 0, ( |3.5| ) is just 
O(o) = d; this was immediate from the definition of ar n ) . 

Let Llj(^) be the orthogonal projection onto the Aj -eigenspace. The IT satisfy the 
conditions 

n? = it, 
n,n fc = o (i^k), 

s 



■ V- 



=1 



tryLIj = aV (3.6) 

In contrast to the eigenvalues, the projections depend on the direction £/|£| of £, rather 
than on the magnitude |£|. In other words, they are O-homogeneous in £. Furthermore, 
they are polynomial in £/|£|: 

2p 
n=0 

for some p, where 

1 

IIP 

Here the n»( n ) are some End (V)-valued symmetric n-tensors that do not depend on £. 

There is, however, quite a bit of ambiguity in the definition (|3.7| ) of the homogeneous 
polynomials Pi( n ){£,) , since multiplication of an n-homogeneous polynomial q(£) by |£| 2 
produces an (n + 2)-homogeneous polynomial without changing the associated 0- 
homogeneous function g(£/|£|). We can remove the ambiguity by requiring that Pi( n ) have 
no |£| 2 factor. This is equivalent to requiring that Pj( n )(£) be a harmonic polynomial in 
£, which in turn is equivalent to requiring that its restriction to the unit ^-sphere is an 
n— -order (End(V)-valued) spherical harmonic. Yet another equivalent formulation is to 
require that IIj( n ) is trace free in all its indices. At any rate, with this requirement, we 
have uniquely defined quantities LL( n ) . 

Note that the explicit formula 

n _ (A - AQ ... (A - ViKj - Ai+Q ...(A — X s ) 
(Aj — Ai) . . . (Aj — Aj_i)(Aj — Aj + i) . . . (Aj — A s ) 

exhibits each projection as a polynomial of degree 2(s — 1). By homogeneity, 

(Ami) - a*o • • • mm) - ^m/m - /*+o • • • mm) - ^ 



n*(0 



(yUi - Ul) . . . (/ii - - fli+l) ■ ■ ■ (fa - Us 
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We can take this expression and expand as a product of homogeneous terms: 

11,(0 = E^^ 1 WW) = E^« (WW • ..a^-^- ^'^ k -\ (3.8) 
k=i k=i 

where are numerical constants depending only on the fij. Since A(£) is 2- homogeneous 
it follows that all homogeneity orders are even, i.e. Pi(2n+i) (0 = n«(2n+i) = 0, and 
therefore, 

v 

n «(0 = E p ^)(0- (3-9) 

n=0 

In turn, by writing A(£/|£|) = Yli=o ^i^iiO we compute powers of A in terms of projec- 
tions: 

s 

^(£/l£l) = E^ n ^)- 

Substituting this into (J3.8|), we get 

s 

fe=l 

In other words, the matrix of coefficients C := (cjjj.) is inverse to the (Vandermonde) 
matrix of powers M. = (ftkj) '■= 

(cm,) := (/%)"* • (3.10) 

Though the number s depends on the particular leading symbol a, there is an upper 
bound for s which depends only on the representation (ip, V) to which the bundle V is 
associated. This is described in detail in ||. In particular, in case ((p, V) is irreducible, the 
algebra generated by restrictions to the unit ^-sphere of equivariant leading symbols of all 
orders is commutative, and thus simultaneously diagonalizable. The resulting projections 
diagonalize equivariant leading symbols of any (not just second) order. The number of 
projections, i.e. the dimension of the algebra just described, may be described in terms 
of representation-theoretic parameters. For reducible (cp, V), similar considerations are 
valid, but the algebra is not commutative. 

Now write the leading symbol in terms of projections 

s s p 

a(o = iei 2 5> n ^) = lerEE^wo- (3.H) 

i=l i=l n=0 
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Since A(£) is 2-homogeneous, it is a sum of End(V)-valued spherical harmonics in £ of 
degrees 2 and 0. Thus the spherical harmonics on the right of all other orders vanish: 

s s 

yilfr-Pi(2n)(Q = ^2^i(2n) = 0, U ^ 0, 1. 

i=l i=l 

Furthermore, it is clear that the decomposition of A(£) into n = 2 and n = contri- 
butions comes upon taking 

Mt) = a'"t& = \t\% + b£'t ll Z v , 

where b^^^u is a second-order spherical harmonic in £; i.e., 6^ is trace free in its two 
indices. This gives 



where 



^ 5 5 



i=l i=l 



4 Symmetric Two- Tensors 

It is instructive at this point to consider two examples: the bundle S 2 of symmetric 
two-tensors, and the subbundle Sq of trace-free symmetric two-tensors. In particular, in 
these examples, one begins to glimpse the differences and relations between the cases of 
reducible and irreducible V. We may compute with tensors valued in either complex or 
real tensor bundles; in fact, most of the following discussion holds in either setting. But 
for the sake of defmiteness, and with a view toward applications, let us assume that all 
our tensors are real. 

First note that a basis of the 0-homogeneous symbols in our class is given by 

(Xl<£>) q/ 3 = g a f3^fi , 

(X 2 y?) a/3 = jjpf "£(0^ » 

(X 3 y?) a/3 = j^Za&^n , 
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Note that in the inner product (ip, ip) = ip^tp^ 11 , we have 

A x — Ai, A 2 — A2, A 3 — A4, A 5 — A5. 

In addition, the traces of these endomorphisms are 

m(m — 1) 
tr 52/52 = i, 

ti s' 2 Xi = m, 

tr 5 2X 2 — m+1, 

tr^A^ = tr 5 2X 4 = tr^Xs = 1. 

The multiplication table of the algebra generated by Iy and the Xj is shown in Table 
1. The arrow indicates that the left factor is to be found in the leftmost column. For 





x l 


x 2 


x 3 


x A 


x 5 


x 1 


mXi 


x. 


x l 


771X4 


x± 


x 2 


x 3 


\{X 2 + X,) 


x 3 


x 5 


x 5 


x 3 


mX 3 


x 5 


x 3 


mX 5 


x 5 




X! 


x. 


x l 


x. 


x 4 


x 5 


x 3 


x 5 


x 3 


x 5 


x 5 



Table 1: Multiplication table 

example, X 3 X 4 = mX 5 . 

The leading symbol of any even-order operator has the form 

-4(0 = l£| 2p ("o Is* +»iX 1 + 2a 2 X 2 + a 3 X 3 + a A X 4 + a 5 X 5 ) , 

for some p and some numerical parameters a,. If A is the leading symbol of a second-order 
operator, then o 5 vanishes. If in addition A is self-adjoint, then 0:3 = 0:4. Thus for a 
second-order self-adjoint leading symbol A, 

A(0 = Kl 2 K I& +a 1 X 1 + 2a 2 X 2 + a 3 (X 3 + X 4 )) . (4.1) 

If we further require that A be positive definite, we get additional inequality constraints 
on the Oj ,; these are discussed below. 
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4.1 Trace-free symmetric two-tensors 

Now consider the bundle Sq of trace-free symmetric two-tensors. The symbol 

P:=I S - -X x 
m 

is the self-adjoint projection onto Sq. Thus we may define a spanning set of symbols Yj 

•2 




on SZ by 



According to Table 1, 



and 



Y = PXiP. 



Y 1 = Y 3 = Y, = 0, 



Y 2 =X 2 --(X 3 + X 4 ) + ^-X l7 
m m z 

Y, = X,--{Xz + X i ) + — 2 X l . 



Clearly both Y 2 and Y5 are self-adjoint. Together with the identity symbol J52 , the 
symbols Y 2 and I5 form a basis of the symbol algebra on Sq . For purposes of explicit 
computation, it is useful to note that 



PX 1 = XiP = PX 4 = X 3 P = 0. (4.2) 



We compute that 



Tfl — 1 

Y 2 Y 5 = Y 5 Y 2 = Yl = F 5 , (4.3) 

m 

2 2 2m K ' 

This points up a very important and useful fact: that the symbol algebra over Sq is 
commutative. In fact, as shown in Jj], this is a general feature of the case of an irreducible 
bundle. In the case of a bundle, like S 2 , that is reducible under its structure group, the 
symbol algebra may be noncommutative. 

The commutativity of the symbol algebra over Sq significantly simplifies the computa- 
tion of the projections. First note that all symbols will be simultaneously diagonalizable, 
so the discussion on this level is independent of the particular symbol A. From ( [4.3|) we 
immediately obtain one projection, namely 

Tfl 

n 3 = rY 5 . (4.5) 

m — 1 
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Since I S 2 , Y 2 , and Y 5 form a basis of the symbol algebra, ( |4.4| ) shows that J52 , Y 2 , and 



Y 2 2 are also a basis. Note that 

3m — 2_ Q m — 1 

V"- 

and 



Yo = Y 2 2 - — — -Y 2 , 

2 2m 2 2m 



If I! := a/52 + bY 2 + cY" 2 , we thus get 



4 _7m 2 -10m + 4 2 (3m-2)(m-l) 

I2 — : — o ^2 a — 5 ^2- 

4m z " 4m 

2 



6c(m — 1) c 2 (3m — 2)(m — 1) 
m Am 2 



n 2 = a 2 I S 2 + {2ab- — ^ - - , ^ n v "" — ^ J> F 2 



f „ 6c(3m-2) c 2 (7m 2 -10m + 4)] 

+ ^ 2ac + o 2 + — i - + — i — ^ i. Y 2 . 

y m Am z J 

There are 8 = 2 3 solutions (a, b, c) of the projection equation n 2 = EL Two of these are 

Aim — 1) Am 



(a 2 , b 2 , c 2 ) := 0, 

\ m — 2 m — 2 

and 

m 2m 2 



(a 3 ,o 3 ,c 3 ) := 0, 



(m — l)(m — 2) (m — l)(m — 2) 
these in fact correspond to fundamental projections EI2 and EI3 . (EI3 is the same as the 
projection found by inspection in Q4.5I) .) 

The remaining solutions correspond to the third fundamental projection Eli := I S 2 — 
EI2 — EE3 , and to 0, I S 2 , 112 + 113, Eli + EI3 , and IIi + EI2 . The fundamental projections 
have particularly simple expressions in terms of Y 2 and Y5: 

n 2 = ^^(-my 2 2 + (m-i)y 2 ) = 2(y2-n) (4.6) 

n 3 = ? ^ -(2Y 2 2 -Y 2 ) = ^-Y 5 (4.7) 

(m — l){m — 2) m — 1 

vn — 2 

IL X = Jo 3 -2Y 2 + -Y 5 (4.8) 

m — 1 

The general self-adjoint second-order symbol A on Sq may be viewed as the compression 
PAP of a second-order symbol i on 5 2 , which, in view of ( |4.1[ ), is 



4>(£) = |e| 2 («o^+2a 2 F 2 
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= |£| 2 (A*iIIi + ^112 + ^3113) 

= |e| 2 {^ 1 / 5 2 + 2(/i 2 -/i 1 )r 2 

^—j [(m - 2)/i 1 - 2(m - l)/i 2 + m/i 3 ] Y" 5 j. 



+ 

m 

From this we obtain the eigenvalues in terms of the parameters a and o 2 : 

2(m- 1) 

/ii = a , A* 2 = « + a 2 , fj, 3 = a -\ a 2 ■ 

m 

Thus, self-adjoint positive definite second-order symbols are in one-to-one correspondence 
with choices of (cto , a 2 ) for which 

a , 0:0 + 0:2, ma® + 2(m — l)o 2 > 0. 

Note that for m > 2, the first and third conditions together imply the second one. 

Our projections have the following interpretation. After choosing the distinguished 
direction £, one can distinguish three subspaces of the trace-free symmetric two-tensors. 
First, there are tensors in the direction of the trace-free part (£ <8> £) of £ <g> £. This is 
clearly the range of II3 , since Y 5 ip is a scalar multiple of (£ <E> £)o- Next, there is the 
subspace consisting of tensors £ V (, where ( _L £. This is the range of n 2 , since for 

C/3 ■= j|j£V/JA ~ ]Jj3 O^^'Va/. ' 

we have 

(y 2 - r 5 V = e v c, C A 6 = 0. 

The remaining subspace is the range of Hi ; it may be described as the space generated 
by tensors (( V 9) , where ( _L £ _L 6*. The traces of the projections are thus 

(m + l)(m-2) „ „ 

tr v n 1 = ^ ^ ^, tr v n 2 = m-l, tr v n 3 = i. 

4.2 Symmetric two-tensors (not necessarily trace- free) 

To leave the irreducible setting, consider the bundle S 2 of symmetric two-tensors (unre- 
stricted as to trace); this is equivariantly isomorphic to the direct sum of the trace-free 
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symmetric two-tensors and the scalars. In fact, the direct sum decomposition is imple- 
mented by the projections P and J52 — P = {l/m)X\ . By noting that X\ = m(I S 2 — P) 
and using Table 1, we obtain in addition to 











m 


^ m 


= ± 








m 


^ m 


= x 3 


- 

m 


= x 4 





1 > 



(J v -P)X 2 (/ v -P) = -(/ v -P), 

m 

(J v - P)X 3 (/ V -P) = (Iv - P)X A (I V -P) = (J v - P). (4.9) 
Let III, n 2 and II3 be the fundamental projections defined by (|4.6| ). Denoting 



T = - r L= (x 3 - -X,) , T* = -=L= (x 4 - -X x 
\Jm - 1 V m J sjm — 1 \ m 



and using the equations (|4.9|), we obtain the decomposition 

A{i) := |£| 2 [/nIIi + /i 2 n 2 + /i 3 n 3 + K (T + T*) + g (/ 52 - P)] , 



where k and g are real constants defined by 



2 



k = \/m — 1 ( — a 2 + «3 
\m 

2 

g = a + mai H a 2 + 2a 3 . 

m 

Note the useful relations 

tp = tit = Pn 2 = TTI3 = 0, 
pt* = n x T* = n 2 T* = n 3 T* = o. 

T 2 = (T*) 2 = 
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T*T = I V -P, TT* = n 3 . 

n 3 T = T, T*n 3 = T*. (4.10) 

Observe also the obvious non-commutativity; the projections will now depend on the 
particular leading symbol we are studying. We see that for our new symbol on S 2 , the 
first two projections IT , Il 2 computed above are untouched, and there are two more 
projections Z 3 , Z± onto one- dimensional subspaces. IT and n 2 are independent of the 
particular symbol we are diagonalizing, while and Z 3 and Z 4 depend on it. Z 3 and Z 4 
take the form 

Z = aU 3 + 6(T + T*) + c{I v - P) 

By ( fllO|) , we have 

Z 2 = (a 2 + 6 2 )n 3 + b{a + c)(T + T*) + (b 2 + c 2 )(/ 52 - P). 
The projection equation Z 2 = Z gives 

(2a - l) 2 + (26) 2 = (2c - l) 2 + (26) 2 = 1, b(a - c + 1) = 0. 
Aside from the solutions 

(a,b,c) = (0,0,0), (1,0,1), (4.11) 

all solutions have the form 

a — -(1 + cos 9), 6 = -sin6 l , c—-(l — cos 9), 

where 9 is an arbitrary real parameter; conversely, all choices of 9 give solutions. The 
solutions ( |4.11| ) give 0- and 2-dimensional projections, so can be discarded from the present 
point of view, where we seek complementary one-dimensional projections. If we define 

(a 3 ,6 3 ,c 3 ): = Q(l + cos0), ^sin0, ^(1 -cos(9) 

(a 4 ,6 4 ,c 4 ): = Q(l-cos6»), -isin6», ^(1 + cosi 

we get a set of complementary projections Z 3 and Z 4 ; that is, we have Z 3 Z^ = Z^Z 3 = 0. 

We still need to find a value of 9 adapted to our given symbol A. Denote by z/ 3 and 
v± the eigenvalues of A in the ranges of Z 3 and Z 4 respectively: 

/i 3 IT 3 + k (T + T*) + q(I v -P) = u 3 Z 3 + u 4 Z 4 . 
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We find that 

(1/3 + ^4) + (^3 - ^4) cos6» = 2/i 3 , 
(1/3 — 1/4) sin^ = 2k, 
(z/ 3 + 1/4) - (u 3 - z/ 4 ) cos# = 2q. 
From these equations we first determine the eigenvalues 

Z/ 3 = p + LJ, V 4 = p - LU, 

where 

1 m 

p = -{fj, 3 + q) = a + —an + a 2 + o 3 , 



2 2 , / /':< - V ' 2 



1 2 

[— m 2 a;i + 2(m — 2)a 2 — 2m« 3 ] + 4(m — l)(2a; 2 + ma 3 ) z 



4m 2 

The positivity of the leading symbol is translated now into the condition uo 2 < p 2 , or 
k 2 < qp 3 . So, in addition to a > and oo + 02 > we have 

(m — l)(2o;2 + rnas) 2 < (ma + m 2 ai + 2a 2 + 2ma 3 ) [ma + 2(m — 1)0:2] • 
The parameter # is now determined by 

cos# = — — - = [— m 2 ai + 2(m — 2)a 2 — 2ma 3 l , 

2uj 2mu L ~ J 

smt) = — = (2a 2 + mas) . 

to mcu 

Therefore, the projections Z 3 and Z4 have the form 
and indeed depend on the symbol A. 
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5 The resolvent and the heat kernel 

Let F be a self-adjoint second-order partial differential operator on a compact manifold 
with positive definite leading symbol. (In particular, F is elliptic.) Then F has discrete 
real eigenvalue spectrum which is bounded below by some (possibly negative) real number 
c. If A is a complex number with Re A < c, then the resolvent (F — A/) -1 is well defined: 
if {(Xj ,<fj)} is a spectral resolution, with the <fj forming a complete orthonormal set in 
L 2 (V), then 

(F-A/)-V, = (A,-A)-V,- 

The resolvent is a bounded operator on L 2 (V), and in fact is a compact operator by the 
Rellich Lemma, since it carries L 2 (V) to the Sobolev section space L\(y) continuously. 
The resolvent kernel is a section of the external tensor product of the vector bundles V 
and V* over the product manifold M x M, and satisfies the equation 

(F-XI)G(X\x,y) = 5(x,y) 

where S(x, y) is the Dirac distribution (which in turn is the kernel function of the identity 
operator). Here and below, all differential operators will act in the first (x, as opposed 
to y) argument of any kernel functions to which they are applied. The resolvent is well- 
defined as long as the null space of F — XI vanishes; i.e., as long as A is not one of the 
Xj . As a consequence of self-adjointness, we have 

G*(\\x,y) = G(\\y,x). 

Similarly, for t > the heat operator U(t) = exp(-tF) : L 2 (M,V) -> L 2 (M, V) is 
well defined. U (t) is a smoothing operator; that is, it carries L 2 sections to sections in 
(IkeitLl = C°°. The kernel function of this operator, called the heat kernel, satisfies the 
equation 

(dt + F)U(t\x,y) = 

with the initial condition 

U(0 + \x,y)=5(x,y), 

and the self-adjointness condition 

U\t\x,y) = U(t\y,x). 
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As is well known [0, the heat kernel and the resolvent kernel are related by the 
Laplace transform: 

oo 

G(X) = J dte tx U(t), 
o 

c+ioo 



U(t) = — [ d\e- tx G(\). 
2m J 



It is also well known [|T^| that the heat kernel U(t\x, y) is a smooth function near the 
diagonal {x = y} of M x M, with the diagonal values integrating to the functional trace: 

Tr L 2 exp(— tF) = / d vol (x) tr v U(t\x, x). 
Jm 

Moreover, there is an asymptotic expansion of the heat kernel as t — > + , 



Tr L 2 exp(-tF) ~ (Ant)- m / 2 ]T ^A k , 



k\ 

k>0 



and a corresponding expansion of the resolvent as A — > — oo: 

Tr L *d n x G(\) ~ (4vr)- m / 2 ^ ^pHC* - m)/2 + n + 1](— A)™/ 2- * -71-1 ^, 

fc>0 

for n > m/2. Here A k are the global, or integrated heat coefficients, sometimes called the 
Minakshisundaram-Pleijel coefficients. 

There is additional information in the local asymptotic expansion of the heat kernel 
diagonal, 

U*v(t) := U(t\x, x) ~ (47rt)- m / 2 ■ (5.1) 

fe>0 

The local heat coefficients a k integrate to the global ones: 



Ak = J d vol (x) tr vct/c- (5.2) 

One can, in fact, get access to the local heat coefficients via a functional trace, by taking 
advantage of the principle that a function (or distribution) is determined by its integral 
against an arbitrary test section / e C°°( End (V)): 

Tr i2 /exp(-tF) ~ (Ant)^ 2 ^ ^Mf, F )- 

k>0 
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We have 



A k = A k (l,F), A k (f,F) = J dvo\(x)f ti v a k . 



M 

The local heat coefficients a k have been calculated for Laplace type operators up to 04 
0. For non-Laplace type operators, some of them are known only in the very specific 
case of differential form bundles |J. 

We shall calculate below the coefficients A and A\ for NLT operators in terms of the 
projections introduced in the previous sections. Our tactic will be to construct an approx- 
imation to the heat kernel U(t\x,y), that is a parametrix. The important information in 
the parametrix for small t is carried by its values near the diagonal, since the heat kernel 
vanishes to order 00 off the diagonal as t — > + . Since the heat and resolvent kernels are 
related by the Laplace transform, this is equivalent to studying an approximation to the 
resolvent kernel G(X\x,y) near the diagonal for large negative Re A. 

Let us stress here that our purpose is not to provide a rigorous construction of the 



resolvent with estimates; for this we rely on the standard references |[L2| . Rather, given 
that the existence of resolvent and heat parametrices is known, our aim is to compute var- 
ious aspects of it; in particular, information on leading order terms sufficient to determine 
some of the heat kernel coefficients A k . 

We shall employ the standard scaling device for the resolvent G(X\x, y) and heat kernel 
U(t\x, y) when x — > y, A —>■ —00 and t — > 0. This means that one introduces a small 
expansion parameter e reflecting the fact that the points x and y are close to each other, 
the parameter t is small, and the parameter Re A is negative and large. This can be 
done by fixing a point x', choosing normal coordinates at x' (with g^ u (x') = 5^ v ), scaling 
according to 

x — > x £ = x 1 + e(x — x'), y — > y £ = x' + e(y — x'), 

t — > t e — £ t, A — > \ E — £ A, 

and expanding in an asymptotic series in e. If one uses a local Fourier transform, then 
the corresponding momenta £ G T*M are large and scale according to 



This construction is standard | 12| . This procedure can be done also in a completely 
covariant way 0, ||] . 
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In the case of Laplace type operators, the most convenient form of the off-diagonal 
asymptotics as t — > 0, among many equivalent forms, is [Q. [3J 

U(t\x, y) ~ (Anty m / 2 exp (-£) A 1/2 £ ^&*(*, y), (5.3) 

fe>0 

where cr = a(x,y) = r 2 (x,y)/2 is half the geodesic distance between x and y, and A = 
A(x,y) = \g\^ 1 ^ 2 (x)\g\^ 1 ^ 2 (y) det (—d^d^,a(x,y)) is the corresponding Van Vleck-Morette 
determinant. The functions bk{x,y) are called the off-diagonal heat coefficients. These 
coefficients satisfy certain differential recursion relations. Expanding each coefficient in 
a covariant Taylor series near the diagonal, one gets a recursively solvable system of 
algebraic equations on the Taylor coefficients @ . The diagonal values give the local heat 
kernel coefficients a^(x) = bk(x,x). However, in the general case of non-Laplace type 
operators, it is very difficult to follow this approach, since the Ansatz for the off-diagonal 
heat kernel asymptotics ( |5.3| ) does not apply; the correct Ansatz would be much more 
complicated. For this reason, we employ the approach of pseudo-differential operators (or, 
roughly speaking, local Fourier transforms). An alternative approach, which generalizes 
the Ansatz fl5T3]), is developed in Section |T0|. 

6 Gaussian integrals 

The Gaussian average of a function /(£) on R m is 

</>= /"^-' s|2 /(fl. 



71 



R m 

The Gaussian average of an exponential function gives the generating function 
h(x) = j exp (-|£| 2 + z£ • x) = (exp (i£ • x)) = exp ) ■ 

R m 

Expansion of Iq(x) in a power series in x generates the Gaussian averages of polynomials 
< 1 >= 1, < ^ >= 0, < >= -g^, 

< SMI ■ ■ S/i2n+l > = 0j 

_ (2tQ! 
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Furthermore, using the relation 

and analytic continuation, one can obtain the more general formulas 



(6.1) 



S>1 £, 



M2n 



r(m/2 + n-p) (2n) 



0, 



|^|2p / " r(m/2 + n) 2 2 »n!^ 1 " a ' ' ' 9 ^~^ 
for any p with Rep < n + m/2. This means, in particular, that 



S/il ' ' ' s 7 , 



/<2n 



lsT p 

SMl ' ' ' S>2n 



r( m / 2 + n-p) ( w-^r, 
r(m/2) \ |e| 2 " 



lei 



2n 



T(m/2) (2n)! 
r(m/2 + n) 2 2 ™n! 



5 I (M1M2 ' ' ' fl , M2n-lM2n) 



(6.2) 



(6.3) 



Note that if a function / depends only on £/|£|, that is, if it is homogeneous of order 
0, then the average introduced above is a constant multiple of the average over the unit 
(m — l)-sphere S™ 1 ^ 1 . 

We will also need to compute Fourier integrals of the form 



'Ml---M2r. 



(x) 



R'-' 



S^Ml ' ' ' s"M2n I •/- \ 

HP exp^-x) 



(6.4) 



The trace of the symmetric 2n-form 7 2 n over an Y two indices is 2(ra — l)-form -?2n-2 : 



n ^2n-l/J.2n J —J 

i/ - I Ml'"M2n-2M2n-lM2n Ml'""M2n-2" 



Therefore, the total trace of the symmetric form 7 2 n is -^o : 



tr fl / 2 n(^) = Io(x) = exp ^— 



By rescaling £ — > V^C one can show that this integral satisfies the equation 



t m/2 d n 



, -m/2 j I 



^mi • ' • d ^ n exp 



\x\ 



At 
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Taking into account the obvious asymptotic condition lini£_ >0O [t m / 2 I n (x / \ft)\ = 0, by 
multiple integration we get 

By changing the variable, s = 1/u, we can rewrite this in the form 

W-^M = T^^y f duu m l 2 - n -\l -ur-%, ■ .-^exp (-^«) . (6.5) 



o 

If m is even this can be computed in elementary functions. 

In the case n < m/2, one can interchange the order of the integration and differentia- 
tion. Then by using the formula 

1 du u a -\l - u) M e™ = T{a) ^~ a) iFa (a; b; z) , 

r(6) 

(with Re b > Re a > 0) where 

k 







^ r(a + A:)r(fe) 
1 F 1 (a-b-z) = ^— — — — z 



[.. _ T(a)T(b + k)k\ 
is the confluent hypergeometric function, we obtain 

\x\ 



where 

~ , , lYm/2 — n) tm m \ v-^ Tim/I — n + k) h 

*"W = r( , n /2) lFl (t " n; T ; "7 = g fR2TW 2 

(with n < m/2). 

If n > m/2, then this formula cannot be applied directly. However, it is still valid if 
one does analytic continuation in the dimension m. The physical (integer) value of the 
dimension should be put after computing the derivatives. Alternatively, one can do the 
differentiation in ( |6.5|) explicitly before the integration. This is equivalent to subtracting 
the first n terms of the Taylor series of the exponential. This subtraction is clearly 
harmless, since the 2n-th derivatives of these terms vanish. However, this makes the 
integral finite, and justifies the interchange of the differentiation and integration. The 
result applies for any n. In this way one obtains 

= (-l) n ^ 1 ---^ 2 A('-M!y (6.6) 
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where is obtained from the function & n {z) by substracting the first n terms of the 

power series: 

. , r(m/2 — n) „ (m m \ v-4 Y(m/2 — n + k) , 
A r(m/2 - n + fc) fc 

^ r(m/2 + A;)A;! ' 



k=n 



7 Leading order off-diagonal heat kernel asymptotics 

In general, the leading order resolvent and the heat kernel are determined by the leading 
symbol of the operator F: 



where £ G T*M, £ • (x — y) = — V^)i an d ^£ is Lebesgue measure on R m . Here and 

everywhere below all integrals over £ below will be over the whole R m . 



Writing the leading symbol in terms of the projections from ( 3.11 ), it is not difficult 
to obtain 

< _ e *(*-i/). "<(o 



The problem is now to compute these integrals. If we are only interested in traces, 
then by Q3.6|) , we can easily compute 

tr v £/o(*|z,y) = Ed i (47rt /Ui )- m / 2 exp (-^-^ 
i=i \ ^ 

That is, the trace of the leading order term in the heat kernel asymptotics is a weighted 

linear combination of the scalar leading order heat kernel term with scaled times, t [lit. 

Similarly, we can relate the trace of the resolvent at leading order to resolvents of Laplace 

type operators: 

tr vG (A|x, V) = J2 di{2itm)- ml2 (^T^2 ) ( ^ K {™-2)/2 (J^\ x ~ v\\ > 



Ivan Avramidi and Thomas Branson: Non-Laplace type operators 



30 



where K p (z) is the modified Bessel function. 



Applying the methods described in the previous subsection we can compute the heat 
kernel and resolvent before taking the trace: 



U (t\x,y) = Y J ^iY m ' 2 (exp 



i=l 



VtJTi 



n*(0 > • 



Using the decomposition of the projections into spherical harmonics ( |3.9| ), we obtain 
U (t\x,y) = Yl (47r^)- m/2 q fc a^^---a^- 3 ^- 2 )(-^) fe_1 

Ki,k<s 



\x — y\ 

x d m ■■ -dun-^k-i ( ^— , • 



(4vrt^)- m / 2 Q fe (t/i i Fo) fc - 1 <l> fc _ 1 



KiKi 



\x-y\ 
4tfa 



(7.1) 



G (X\x,y) 



,-m/2 



k-l 



l<i,k<s 

(•DO 



X 



(ite tA t -m/2+fe-l $fc _ i 



F _ 1/1 



'o v 
where F = a^d^dy and the functions $ n (z) are defined in previous section; they are 
given by ( |6.7| ). We would like to stress that these formulas can be presented locally in a 
"covariantized" form. This is effectively achieved by replacing \x — y\ 2 by 2o~(x,y) and 
(x — yY by V^a(x,y) and adding a factor A 1 / 2 (x,y) (for details, see 0); the objects a 
and A being defined after equation (|5.3|). 



8 The heat kernel coefficient ao 

The leading order term of the diagonal heat kernel asymptotics can now be written as 

Ut S (t) := U (t\x,x) = (47rt)- m / 2 ^/i- m/2 (IT) . 
This gives the coefficient ao : 



i=i 



(8.1) 
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The Gaussian average is computed by using (|6. 3[) : 

r(m/2)(2fc-2)! 



<n i >-n 2(0) - l.^ r(m/2 + A; _ i)22fc _ 2(fc _ i)! tr 9 V a, (8.2) 

where V n is the symmetrized tensor power of symmetric form, tr g is the total trace of a 
symmetric form defined in (3^) and the constants are defined by ( gjjj) . 
From this, we have the formula 

Ed% 

8=1 A 4 ! 

for the trace of the heat coefficient. 

Note that if our bundle V is associated to an irreducible representation of Spin(m), the 
averages of the projections, being Spin(m)-invariant endomorphisms, are proportional to 
the identity, by Schur's Lemma. The exact proportionality constant may be determined 
by taking the trace. We have: 

di „ 1 di _ 



TT d% 1 \ "V 

<il>=-/ v , «o = ^2^ 



d u d ^ u ™'* v - 

i=l H'l 

These formulas point up a new feature of non-Laplace type operators; one which 
complicates life somewhat. Whereas the dimension dependence of the heat coefficients 
of Laplace type operators is isolated in the overall factor of (47r) _m / 2 , the dimension 
dependence for NLT operators is more complicated. 



9 The heat kernel coefficient A\ 

As we have seen, even the computation of the leading order heat kernel requires significant 
effort in the case of non-Laplace principal part. This indicates that the calculation of 
higher-order coefficients will be a challenging task. In this paper we will compute the 
coefficient A\. Since A\ is the integral of tryai by ( |5.2|) , it suffices to compute the local 
coefficient a\ modulo trace-free endomorphisms. By elementary invariant theory, a\ has 
the form 

ai = Yl ^(fcV™^ 1 )''^ + H X R + H» V R W + H^R^p (9.1) 

where Q/w, P(k), an d H% are End (V)-valued tensors, and R, R^ u , and R^ ua/ 3 are the 
scalar, Ricci, and Riemann curvatures. 
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For the trace we have a similar formula, 

tryai = tv v (h q) + h x R + h^R^ + h^ uaf3 R pua p 

where hi = tryfZ, (i = 1,2,3) are some tensors. By invariant theory, we may conclude 
that these tensors have the following form: 



ho - ^ k Q(k)i 



fc>0 



h = d, /if = C2ST, h^ aP = csg^gr" + c^gT + c^Ttf*, 

where q are some constants. 

Note that for irreducible representations the endomorphisms Qo, Pq, ho and H\ are 
proportional to the identity endomorphism; in particular, h = Coiy, Hi = {c\/d)Iy. 

Taking this into account, we obtain 

tr v ai = tr v(h q) + (3R, 

where 

{3 = ci + c 2 + c 3 - c 4 . 

Thus it suffices to compute the constant (3 and the endomorphism ho- 

To compute ho, we can employ a by now standard variational principle. If the operator 
F(e) depends (in a suitably estimable way) on a parameter e, then 

^ Tr L 2 exp(-tF) = -t Tr L2 (<9 £ F) exp(-tF). 
If our variation comes from scaling q, 

q -> eg, 

then (d e F) = q. Expanding both sides in powers of t and comparing coefficients of like 
powers, we obtain 

d £ Ai= I d vol (x) tr y(a g), 
A/ 

so that 

Ai= / dvo\ (x) [tr v (a q) + f3R] , 
and, therefore, h = a , where a is given by ( |8.1| ) and 
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The computation of the constant (3 is considerably more difficult. It is clear the (3 is 
universal to our class of operators; thus we may compute it for any particular operator 
and manifold, or class of such. Accordingly, note that (3 is equal to the coefficient a\ when 
q = and R — l: 

(3 = trydi 

q=0,R=l 

The following considerations will be completely local. Fix a point x', and compute 
in normal coordinates centered at x', with 9^(0) = 5^ v . Furthermore, impose by gauge 
transformation the Fock-Schwinger gauge for the connection 1-form A. We then have 

[g^{x) - = 0, A^x)^ = 0. 

Further, we can expand all quantities in Taylor series about x' and restrict our attention 
to terms linear in the curvature. This gives 

9fJ,is( x ) &iw T^R/iai/fS^ ^ + . . . , 

det g^{x) = 1 - ^R af3 x a x p + 

Here and below, the dots denote higher-order terms in the curvature. TZ^ is, of course, 
the curvature of the bundle V, given by TZ^ = R a p^T 13 a . The Taylor expansion of the 
leading symbol section is determined by the equation 

V M a a/3 = 0. 

From this, we get 

cr(x) = aT + \a x ^R u \ xp x a x p + .... 

Here and below, we denote a^ u (0) simply by a^ u . By the above, the potential term q will 
only enter the calculation through q(0), which we denote simply by q: 

q(x) = q + 

Since the constant (3 is universal, we are free to compute in the case of a constant 
curvature metric, i.e. 



R 



R 



pua/3 



m\m 



'—^{9 m gpv - gppgav), (9.2) 
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R 

R^u = —g^u, R= const . 

771 

Now let us take the total symbol of our operator in normal coordinates, a(F\x,£), and 
expand in a Taylor series: 



where 



a(F\x, = a L {F\0, f) + a{F x \x, £) + . . . , 

a L {F\Q,t) = ar£ l £ v = A{t), 
a(F 1 \x,0 = -F/^ + .YV^?, 

X^ap = —-a x< - p R u \ a \x\i3), 

and [A, B] + = AB + BA denotes the anticommutator. 

There are many equivalent ways of constructing the heat kernel asympotics on the 
diagonal locally Using the Volterra series 

exp(-tF) = exp(-tFo) - t I drexp[-t(l - r)F ]Fi exp[-trF ] + . . . . 

Jo 

and the formula for the heat kernel on the diagonal 



x=0 



we get 



f/diag^ = J | e _ U (£) _ t £ dre -t(l-T)A®P ie -tTA® + y 

where 

F 1 = a{F x \id^ = X^apd^dl^ - + g. 

Finally, by scaling the integration variable £ — > t~ 1//2 £, we obtain the standard asymp- 
totic expansion of the heat kernel on the diagonal (|5.1f) , with the coefficients 



a = / — iKe' m , 



ai 



f * f'dre-^-^F.e-^. 
J Jo 
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In particular, we recover the formula for ao derived in the previous section. 

To integrate by parts in this formula, we need to know how to differentiate the expo- 
nential e A . This can be done via the Duhamel formula 

Jo 

where 

The contraction of this with £ leads to a much simpler formula: 

Z a d%e- TA ® = -2TA(0e- TA{ °. 
The exponential itself is computed by using the projections II, : 



s 



e -^(0 = £ e -r W i€i a n.(£). 
i=i 

Note that when computing ^-derivatives one can, using integration by parts, act in either 
direction. This trick can be used to avoid having to compute the second derivative of the 
exponential e A . First, we rewrite F 1 in the form 

Fi = d^X aP dl + d%L a + L a dl + Q, 

where 

XapiO = X^ u a p^ u = --a x(M R v \ a \ X \^^u, 

L a (0 = If^, 

11 1 

L^ a = —-a^ x R\ a + —a^R^xau + -R p CTav [T r7 p , a^ v } + . 

Q = q- \a^R, v . 
o 

Now, integrating by parts over £ (and changing r to 1 — r in some places) we get 

x J a e~ SlA X a pe~ {T - S2)A J^e~ S2A 

+2 I dr [ ds \e-^ A J a e- sA L a e- (1 - T ^ A 
Jo Jo 

-e'^ A L n e~ {T - s ^ A J a e- sA 



+ £ dTe-^-^Qe-^. 



(9.3) 
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By separating the curvature factors in ( |9.3| ), one can compute each of the tensors Hi 
entering ( |9.1| ). Note that for a Laplace type operator (the case a^ v = g^ u I) we have 
A = J" = e\ X aP = -(l/3)i?^^^, L a = -(1/6)R^, Q = q- (1/6)RI. 
Therefore the first two terms (with X and L) vanish in the Laplace type case, and we get 
the well known result 

ai = q - -RIv- 
6 

Note that the tensors Hi depend only on the leading symbol, i.e. on a^ u . In principle, 
it is possible to compute them explicitly by using the representation of the e in terms of 
projections and Gaussian averages. 

We shall not do this explicitly, but rather compute only the trace of a\. The number of 
projections involved in this calculation is less by one. Computing in the case of constant 
curvature (|9.2|) , and attaching a tilde to X, L and Q in this case, we have 

R 



3m (m 

Ylm{m 



L a (0 = —— — 7— Tr {ata - (3m - 2) J a + 6[T M , J"] + } , (9.4) 

\2m\m — 1) 



n R - 
Q = g- r— a, 

6m 

where 

a = g^a^. 

Now taking the trace and changing s to r — s in the second integral, we obtain 
tr v oi = tr v [ -^W (-4 / dr I ds x [ ds 2 e~ {T - S2)A 



71 



m/2 



JO JO 



+2 [ dr [ ds \e- {1 ~ s)A J a e~ sA - e~ sA J a e- {1 - s)A ] L a 
Jo Jo L 

+e- A Q\. (9.5) 



In the second multiple, the s-integration may be accomplished explicitly, to give 



7T 



1 rl-r 



tryoi = tr v / <|-4 / dr / tZs a / d S2 e- (r - S2)A 



^0 ^0 



X 



jPe-( 1 - T - 8l+ ' a)A J a e-' lA X. 



a/3 
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1 

+2 / drr 
'o 



+e Q 



-TAja e -{X-T)A _ e -(l~r)Aja e -Tj 



(9.6) 



Now let us consider the different contributions separately. 
Q contribution. 
The Q contribution has the form 



tr- 



J 1=1 J 



By changing the integration variable £ to (//j) 1//2 £, we obtain a Gaussian average: 
trv^A*7 W/2 < n i > f?~ 7^~ SjR ) = trva ° y~ 6m^ / ' 

8=1 \ / \ / 

L contribution. 

This contribution has the form 



2tr- 



Using (|9.4j ), we get 



7T 



m/2 



dr r 



-tA ja e -(l-r)A _ e ~(l-r)^ J" e ~^ 



(9.7) 



tr y / — / (h t 



6m(m — 1 



7T 



m/2 







-Ty4 



ja e -(l-r)A _ e -(l~r)Aja e -T/ 



(a£ a - (3m -2)J a + 6[T [va] , J u } + ) 



(9.8) 



Noting that J Q £ Q = A, we see that the term in L a proportional to £ Q does not contribute. 
The term proportional to J a gives 



R (3™ -2) /' < 

6m(m — 1) 



7T 



m/2 



ja e ~(l-^)^ _ e -(l-T)ija e -ry 



rv • 



This is computed by again introducing the projections H : 



n (3m -2) v 

,,fc<s ' 



6m(m — 1) ^— 

l<i,A 



drr 



[MiT+Mfc (l-r)]|$| 2 



_ e -[MfcT+A tl (l-r)]|?| 2 



tr yllj J Q rifcJ Q 
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By changing the integration variable £ to [//jT + /ifc(l — r)] _1 / 2 £, we may express this 
terms of a Gaussian average: 

R-}—, It ^jfe tr v (njj a n fc J a ) , 

6m m — 1 

l<i,fc<s 

where 

«* = T dr r {[^r + Ml " r)]-( m+2 )/ 2 - [// fe r + ^(1 - r )]^ m+2 )/ 2 } . 
./o 

Clearly /c^ = — Kk% ; hence, for i = k 

K>ii = 0. 

For % ^ fc we compute 

TV /o 1 \ ( —m/2+1 —m/2+1 

r(m/2-l) m / _ m/2 _ m/2 \ ^ 



r(m/2 + - [ 2 v» ™ v ^ - ^ 

Similarly, the contribution of the term with T is computed to be 

e -rAja e -(l-T)A 



tr v [ — — [ (It t 



m(m — 1) J 7r m / 2 7 
-e-( 1 -^J Q e-^][T M ,r] + 

' 1 o 



m(m 

l<i,K<S 



R 

^trv 

l<i,fc<s 

Thus the total contribution of L is 



mym — 1 



^ 6m J_ n tr v £ ^.{(3m-2)(n i j-n fc j Q ) 

' ' Ki,k<s ^ 



-6 (( riu .ru k ) + (n, j Q n fct r» t {i 



X contribution. 
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This term has the form 



X 



jP e -(l-T-s 1+ s 2 )Aja e - Sl A _ 



We consider first the term proportional to A. Because A commutes with the exponent, 
we have 

*M^T) "v / ^ jf * C dSl [ Ae ' {T ~ 

xJ a e-^- T - Sl+S2)A J a . (9.10) 
Introducing the projections again, we get 



Ki,k<s 



xe -[(r- S2 + Sl )M l + (l-r- S i+ S2 K]|C| 2 /i .|^|2 n .ja nfe J a> 

Scaling, £ — > [(r — s 2 + SiVi + (1 — r — S! + s 2 )/U fe ]~ 1/ ' 2 £, we obtain 

3m m — 1 x ' 

V 7 l<j,fc<S 

where 

Pik = dr ds 1 ds 2 [(r - s 2 + Si)^ + (1 - r - si + s 2 )p 
Jo Jo Jo 

By changing r to 1 — r and switching the roles of s± and s 2 , we see that 

Pik = Pki- 

For i = k, we easily obtain 



(m+ 4)/2_ 



1 -m/2-2 



For i ^ we compute that 



r(m/2 - 1) J / _ m/2 _ m/2 \ 

-m/2+1 -m/2+1 n 

+2 fh 1 (911) 
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Furthermore, the term proportional to £ in the X contribution is 



JO JO 



e -(r-s 2 )Aj(5 e -(l-T-s 1+ s 2 )Aja e - Sl Aj (a ^ 



Remembering that J a ^ a = A and changing the variables of integration, we have 

4 



-R- 



3m(m 

Introducing the projections IT , this becomes 



tr v C dr-r 2 Ae- TA J a e-^ A J a 

1) J ir m ' 2 Jo 2 



l<i,k<s ' 

Finally, scaling £ -> [r^ + (1 - t)^]- 1 / 2 ^, we get 

4 

— tr v 

l<i,k<s 

where 



^ 3m J_n trv £ 7^<l£l 2 n^n fc j a ) 



For i = fc, 



7« = 



1 -m/2-2 



For i j£ k, we compute that 



r(m/2-i) r -m/2-i o^i" 172 

7ifc = ^7 — /0 ■ oW r\-m[m -2)^ -4(m-2)- 



8r(m/2 + 2)(^- / u fc ) \ ^ ' > 4 -^ 



-m/2+1 -m/2+1 

Thus, the total X contribution is 

4 

— trv 

l<i,fc<s 

where 

&ik Pik 7ifc 1 



i? Q ( A rr tr v V /x i£ r iJfc (|£| 2 IT J a II fe J Q > , 
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Clearly 

an = 0. 



For i 7^ k, 



8r(m/2 + 2)[ (/ii-/i fc ) 3 v ' (/ii-/i fc ) 2 

-m/2 -m/2-1 n 

+4(m-2)-^ + m ( m _ 2 )-^ -I. (9.12) 

Result. 

Summing up the various contributions, we get the main result: 



tr ydi = ^2 m/2 tr y < Ilj > g + (3R, 



where 



1 s 



6m 



1 

L<i,fc<s; i^fc 



,r V /€ i J(3m-2)tr v (n i J a n fe J Q ) 
6m(m — 1) 1 



-6tr v [{m^Uk) + (iv^il^)] t [v 



— ^2 PkViktov (\€\ 2 ILiJ a ILkJ<x) 



3m m 

l<i,k<s;i^k 



with the constants Kjfc and o"jfc given by (|9.9| ) and ( |9.12| ) 



Next we need to compute the Gaussian averages. This can be easily done by usin^ 
the formulas ( |6.2| ). First note that the radial part can be always separated by 



This gives 



< UiJ a u k j >=j< iw a n fc j^ >, 
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where 



< ru k j^ >= f < 3~nJ% >, 
< |ei 2 n,rn fe ^ >= m(m / 2) < >, 



~ \t\ ' ICI 



This allows to simplify the result somewhat: 
= -^E^ m/2 tr v <n 4 >a 



i=i 



^—^ Yl M 3m - 2 ) + 4 ( m + 2 )w<*] tr v (n, J a u k J a 



12(m 

l<i,k<s;i^k 



1 

TV / J ^ifctry 



2(m 

l<i,k<s;i^k 



Tr, 



Note that for irreducible representations the endomorphisms < Ilj >, a — , and 
< UiJ a Uk J a > are proportional to the identity Jy- 

Finally, for the sake of completeness, we list the averages explicitly. By using the 
formulas (lO) we get 



~ ~ x - r(m/2)(2n + 2j-2)! 

1 a k 13 " ^ T(m/2 + n + j- l)2^~\n + j - 1)! 



x 



9(fJ,1^2 ' ' ' 9^2n-3^2n-29v 1 V2 ' ' ' 9u2j-^,U2j-29'y5) 



xa (mM2 ■ ■ ■ a^ n -^ 2n - 2 o? a a VlU2 ■ ■ ■ a v ^-^- 2 a 5) p . 

Denoting by V the symmetric product of symmetric forms and by tr g the total trace of 
a symmetric form, we can rewrite this in the compact form 

~ ~ x - r(m/2)(2n + 2j-2)! 

1 a k p ^ r(m/2 + n + j-l)2 2 ™+^- 2 (n+j-l)! 

Kn,j<s v ' 

xtr 9 [(V n_1 a) V a a (y j ~ l a) V , 

where a a is a vector defined by a a = a^c^. 

This completes the general calculation of heat kernel coefficient tr yc^. More explicit 
formulas may be obtained in particular cases by using explicit formulas for the projections. 
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10 The covariant semi- classical approximation 

In this section we show how the standard semi-classical method can be adapted to compute 
the small-t heat kernel asymptotics of a non-Laplace type operator. In treating the semi- 
classical approximation, we shall follow [19|. 

The object of study is the fundamental solution of the heat equation 



for a non-Laplace type operator 



^d t + F^j U(t\x,x') = 



F = -a^V M V, + q. 



with a singular initial condition 



U(t\x,x') 



5(x, x'). 



;io.i) 



Here e is a small formal parameter. 

Let a^' (x, x') denote the tangent vector to the geodesic connecting the points x and x' 
at the point x', the norm of which is equal to the length of that geodesic. In this section 
we describe a systematic method for constructing the local formal asymptotic solution of 
the heat equation as e — > 0. The initial condition suggests the following Ansatz: 

U(t\x,x') = J(t\x,x') exp (~S{t\x,x')j Q(t\x,x'), 

where S is a scalar function, J is another scalar function defined by 



J(t\x, x') 



1 



det 



1 



2ne 



V^ u/ S(t\x,x') 



and fl(t) has the expansion 

n(t\x,x') ~ J2 £k Mt\x,x') 



y/det g^ix') 



k>0 



in powers of e. The leading asymptotics as e — > require 



S + a^S.^S.y i (f> 
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where S = d t S and S- ifM = V^S. From this it follows that the function S is determined by 
the Hamilton- Jacobi equation 

det v [d t S + d w S ;tl S ;u ] = 0. 
Recalling the eigenvalues of the leading symbol, we have 

s 

WldtS + ^s.^f =o. 
i=i 

This has s different solutions, one for each eigenvalue, determined by 

-^+^^ = 0. (10.2) 

These solutions differ by scaling t — > U = /lit, i.e. S(t\x,x') = So(/j,it\x,x'), where So is 
determined by the equation 

^0+^0^ = 0. 

This is the Hamilton- Jacobi equation for a particle moving in a curved manifold. There is 
a Hamiltonian system that corresponds to each Hamilton- Jacobi equation ( |10.2| ). These 
Hamiltonian systems describe the geodesies parametrized by tj = t/ij. The solution of 
this equation is given by the action along the geodesies connecting the points x and x' 
and parameterized so that x(0) = x', x(t) = x. Thus 

S Jt\x, x') = ^1^1 SMx, x') = ^2^1 

where a(x,x') is half the square of the length of the geodesic. This means that 

Ji(t\x,x') = (47rt/iie)" m A(x,x'), 

where A is the Van Vleck-Morette determinant. As an immediate consequence of the 
Liouville Theorem, this satisfies the equation 

-d t J i + 2V lt (Sl li Ji) = 0, 

or 

(D + a%) J]' 2 = 

where 

D = td t + a^V„. 
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Here and below we denote the derivatives of a simply by adding indices to it, i.e. a h 



V M a, cr Mi , = V^V^o", etc. Recall that 



DA l/2 = i(m - a%)A 1/2 . (10.3) 

Thus the semiclassical approximation as e — > polarizes along different eigenvalues. 
That is, all quantities become dependent on the eigenvalue, and the total solution is the 
superposition of all particular solutions for all eigenvalues. Thus, our final Ansatz is: 

U(t\x,x') = V(47rt^)- m / 2 A 1 / 2 exp (-^—) ^i(t\x,x'), 
Qi(t\x,x') ~ ^t k e k <p [i)k (t\x,x'). 



k>0 



For the function Qi(t), we get a transport equation 

! N i + -L i + M)n i (t)=0, 



where 



eH 2 ' ' et 

1 / I 



Ni = — I a - —a^a^o-v , , 
2jii V 2/Xj J 

U = td t -^ + Ia-VVV^A 1 /* + J- 
2 /ij 2/x. 



M = A-^FA 1 / 2 . 

The initial condition for Qi(t) is determined by the diagonal value of the heat kernel 

Qi{0\x,x) =<Ui>, (10.4) 



where < IT > is defined by ( |8.2p . While for irreducible representations this is proportional 
to the identity matrix, for a general reducible bundle it is not. By using the expansion of 
fij(t), we get recursion relations 

^V*0(t)Q = 0) 

N4 {i)k = -(Li + k — - M^(j) fc _2 k>2. (10.5) 
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Note that Ni is just an endomorphism (i.e., has order zero as a differential operator). 

Now we have a quadratic form a^cx^ov which can be expanded in terms of the same 
projections as before. The role of the covector argument £ of each projection is now played 
by Vcr, and we denote 



■fl2n °JH ' ' ' <T M2n 



(2a) n 

Cika<>"i» ■ ■ ■ a^-3^- 2 )^_^zl , (10.6) 



n=0 

s 

, cr,,,i • • • a 



k=l 



a^o^Oy = 2cr ^ /^P* . 
Recall that V \a^ v = 0; this implies that the tensors Hi( n ) (not Pj) are covariantly constant: 

vn i(n) = o. 

Note that this does not imply that the Pi are covariantly constant, since Vcr is not. Now, 
by using the decomposition of the leading symbol, we observe important properties 

a 

P k Ni = NiP k = — (^ - /j, k )P k , 
^% 

PiNi = N t P t = 0, (/ - Pi)Ni =Ni(I- Pi) = Ni , 

which imply that 

Ni= Yl Y2^i-Vk)Pk- (10.7) 

Q<fe^i<s ^ i 

Next, decompose (fi^k according to the projection Pj : 

<f>(i)k = i/>(i)k + X(i)k, 

where 

1p(i)k = Pi4>(i)k, X(i)k = (I — Pi)4>(i)k- 

Then the recursion ( |10.5| ) takes the form 



NiX(i)o = 0, 



NiX{i)k = ~(Li + k - l)(V>(i) fc -i + ~ M(tfj ii}k _2 + X(i)fc-2) • 

(10.8) 



Ivan Avramidi and Thomas Branson: Non-Laplace type operators 47 



Now, by multiplying this recursion by P n , n^i, and using ( |10.7| ), we obtain 



a 

— 2 (Hi — Vn)PnX(i)l = ~ PnLi1p(i)o , 

-^{Hi - Hn)PnX(i)k = -Pn^{i)h-1 ~ Pn{L>i + k - l)x(i)fc-l 

-P n M(^(i)fc_ 2 + X(i)fc-2) • 

(10.9) 

This recursion determines X(i)k algebraically in terms of V'mjfc-ij V'(j)fe-2, X(i)Jfc-i an d X(i)fc-2- 
In particular, we find that 

X(*)o = 0, 

2 A*? 



X(t)l = ^2 ~~Z P nU^{i)Q- 



l<n=f=i<s 



The recursion does not determine the ip(i)k however. These are determined by another 
differential recursion that is obtained by multiplying ( 10.8| ) by P, , 



PL^(i)o = 0, 

(PiLi + k)ip{i) k = -PiLiX{i)k - PiM(-0(»)fe-i +X(i)k-i), k>l. 
Now let us compute the operator PiLPi entering this recursion. We have 

P^ = td t + -PA-^VV^V.A 1 / 2 

{Pia^a^ + Pa^V^P) - \mP % (10.10) 
where P^ v = V^p. Further, noting that 

Pia^a^Pi = 2/i^P , 

we get 
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Taking into account the equation for the Van Vleck-Morette determinant ( |10.3| ), we obtain 

P^Pi = PiD + K i} 

where 

D = td t + 

is a first-order differential operator, and 

Ki = 1 ^-P l [2a^a^P t;u + (a"" - lug^M P, 

is some endomorphism. Note that <t m = tdx^/dt. Therefore, the operator D is expressed 
in terms of the operator of total differentiation along the geodesies, 

d 

D = t—. 

dt 

The operator D commutes with any function that depends only on "angular" coordinates 
I y/a] in particular, it commutes with the projections: 

DPi = P t D. 

Thus the recursion for the ?/>'s takes the form 

(D + Ki)^ = 0, 

(D + k + Ki)ip(i) k = -PiLiX(i)k — PiM(if>($k-i + X{i)k-i)i k>l. 

(10.11) 

But these are exactly the transport equations along geodesies. They may be integrated 
with the appropriate initial conditions determined by (|10.4j) . In particular, the first 
coefficient ip^o is 

ip(i)o = exp (- J -pK^j < n ; >, 

where the integration is along the geodesic connecting the points x' and x, parametrized 
so that x = x(t) that x(0) = x', x(t) = x. In the standard case, i.e. for a Laplace type 
operator, we have = 0. Therefore, the first coefficient is just 0o — I- 

Thus, the recursion relations determine the asymptotic solution completely. The al- 
gorithm is the following: at each step determine first the X(i)k by ( |10.9|) , and then ipu^ 
by (MM- 
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11 Concluding remarks 

Let us summarize the results of this paper. We have studied in detail a general class of 
non-Laplace type operators, i.e. elliptic second-order partial differential operators acting 
on sections of a tensor-spinor vector bundle over a compact manifold without boundary. 
The only essential assumptions that have been made are: i) the positivity of the leading 
symbol, a^^^ > (in the sense of endomorphisms) for £ ^ 0, and ii) the covariant 
constancy of the tensor a^ v , i.e. Va = 0. We constructed the leading order resolvent 
and the heat kernel and computed the first two coefficients of the heat kernel asymptotic 
expansion explicitly. In the last section we developed an alternative approach for the 
of]- diagonal heat kernel asymptotics for by making use of the general theory of the semi- 
classical approximation. This enabled us to construct a new ansatz for the heat kernel, 
as well as to find a complete set of recursion relations for the coefficients of the off- 
diagonal asymptotic expansion. This generalizes a well known ansatz for the heat kernel 
of Laplace type operators (e.g. see 0). In contrast to the Laplace type case, the off- 
diagonal heat kernel for the non-Laplace type case exhibits some essentially new features, 
notably polarization along the different eigenvalues of the leading symbol. As an explicit 
example of a non-Laplace type operator we considered the most general second-order 
operator acting on the bundle of symmetric two-tensors. We computed the eigenvalues of 
the leading symbol, the multiplicities and the corresponding projections. 
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